Absolute detector calibration using twin beams 
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A method for the determination of absolute quantum detection efficiency is suggested based on the measure- 
ment of photocount statistics of twin beams. The measured histograms of joint signal-idler photocount statistics 
allow to eliminate an additional noise superimposed on an ideal calibration field composed of only photon 
pairs. This makes the method superior above other approaches presently used. Twin beams are described using 
a paired variant of quantum superposition of signal and noise. © 2013 Optical Society of America 
OCIS codes: 270.5570,190.4410,270.5290 



The first suggestions to use photon pairs in absolute 
detector calibration have occurred soon after the exper- 
imental evidence of emission of photons in pairs in the 
process of spontaneous parametric down-conversion had 
been given [1]. The suggested method is based on the 
fact that both photons are created in the nonlinear pro- 
cess simultaneously. Provided that one (signal) photon is 
detected at its (signal) detector, we know for sure in the 
ideal case that the second photon exists in the beam [2] . 
Thus it impinges on a (idler) detector that is calibrated 
with probability one. Many repetitions of the experiment 
then provide the required absolute quantum detection ef- 
ficiency (QDE) [3]. Following this simple scheme, abso- 
lute QDE rji of an idler detector is reached along the for- 
mula rji = (m s m,i)/(m s ), where (m s mj) gives the num- 
ber of coincidence counts at both detectors and (m s ) 
determines the number of signal-detector single counts. 
Weak photon fields having only a small fraction of a pho- 
ton in a detection window on average are needed in this 
approach. The method automatically gives also QDE of 
the signal detector. It can be used for the calibration of 
both analog and photon-counting detectors with the pre- 
cision comparable to other metrology methods exploiting 
thermal and semiconductor detectors [4]. 

The presented method, however, requires as 'the 
probe' a weak field composed only of photon pairs. It 
is unable to cope with any additional noise in the form 
of un-paired single photons superimposed on the paired 
field. A simple analysis shows that the presence of addi- 
tional single-photon noise counts {mi) no i Se modifies the 
formula for QDE rn\ r)i = (m s m.;)/ (m s ) - (m,)^. How- 
ever, if single-photon noise is present in both fields, we 
cannot partition it from the paired part of the probe 
field. The improved formula QDE rji is thus not useful. 
Single-photon noise can only be eliminated as much as 
possible using a careful experimental arrangement. 

Despite this drawback, the method has been gener- 



alized to include also stronger photon-pair fields [5] in 
which both pump-field intensity fluctuations and trans- 
verse correlations of photons in a pair play an important 
role. If the ratio of signal and idler QDEs is known and 
only a paired field is assumed, a simple relation between 
QDE and noise-reduction factor (quantifying sub-shot- 
noise photon correlations) can be revealed [5] . 

As we show in this letter, additional single-photon 
noise in the probe field can be eliminated considering 
photon-number resolving detectors like intensified CCD 
(iCCD) cameras [6,7], array detectors [8] or electron- 
multiplied CCD (EMCCD) cameras [9]. This then al- 
lows to determine both signal and idler QDEs with, in 
principle, no precision limitations. The method is based 
on the measurement of joint signal-idler photocount dis- 
tribution (JPCD) [10]. In detail, the experimental JPCD 
is fitted assuming a special form of the probe field derived 
in the framework of superposition of signal and noise ap- 
plied to paired fields [11,12]. This fitting uses both first 
and second photocount moments and minimizes declina- 
tions from the experimental photocount histogram [13]. 
This approach provides signal and idler QDEs, together 
with parameters of the noisy single-photon signal and 
idler fields giving a detailed detection characterization 
[14]. 

We demonstrate the general method by considering 
the measurement performed by an iCCD camera An- 
dor. Multi-mode twin beams at the wavelength around 
560 nm filtered by spectrally rectangular 14-nm wide 
(FWHM) interference filter have been generated by the 
pulsed third harmonics of a Tksapphire laser tuned at 
840 nm in a non-collinear type-I interaction in a 5-mm 
long BaB204 crystal (for details, see [7]). Whereas the 
signal field has been directly sent to the photocathode 
of the camera, the idler field has been reflected on a 
dielectric mirror (reflectivity 99.2%) first and then im- 
pinged on a different area of the photocathode. This 
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causes asymmetry in the signal (rj s ) and idler (77^ ) QDEs. 
Using pulsed pumping and repeating the measurement 
~ 10 5 times, we have arrived at histogram f(c s ,Ci) giv- 
ing the number of measurements with the registered c s 
signal and Ci idler photocounts. In the experiment, also 
the level of dark counts has to be monitored in order to 
allow for subtracting this additional, but independently 
quantified, noise. 

A suitable choice of a general six-parameter form 
of the joint signal-idler photon-number distribution 
(JPND) together with appropriate values of signal (rj s ) 
and idler (r)i) QDEs lies in the heart of the method. 
However, according to the theory of measurement, only 
values of the first and second moments of measured quan- 
tities are reliably determined after a reasonable number 
of measurement repetitions. Including the five measured 
values of the first and second moments, three free param- 
eters remain in the considered eight-parameter problem. 
Our investigations have shown that the values of remain- 
ing three parameters are successfully revealed requiring 
the best fitting of the theoretical JPND to the experi- 
mental JPCD. 

In the first step, we determine the first ((c s ), (ci)) and 
second ((ci), (c s Ci), (c 2 )) moments of the measured num- 
bers of photocounts. Knowing the first ((d)) and second 
((d 2 )) moments of the dark-count distribution, its con- 
tribution to the measured photocount moments can be 
eliminated and moments of detected photoelectrons can 
be, step by step, found (a = s,i): 



(m a ) = (c a ) - (d a ), 
(ml) = (c? a )-2(m a )(d a )-(dl), 
(m s mi) = (c s Ci) - (m s )(di) - (mi)(d s ) 



(1) 
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On the other hand, the probe photon field 'in front 
of a detector' can be considered as composed of three 
independent parts: pairs, signal noisy photons and idler 
noisy photons. They can be characterized by their num- 
bers M p , M s , and Mi of equally-populated modes and 
mean photon numbers b p , b s , and bi per mode, respec- 
tively. The first and second photon-number moments of 
these fields can be expressed as (a = p, s, i): 

(n a ) = M a b a , ((An a ) 2 ) = M a b a (l + b a ). (2) 

QDEs 7] s and r\i provide the bridge between the 
'theoretical' photon-number moments and experimental 
photoelectron moments. Quantum detection theory [15] 
gives us (a = s,i): 



Va [(n p ) + (n a )] 
((An p ) 2 ) + ((An a ) 2 ) + 

((An p ) 2 ) 



(■m a ), 



f]a 



1<X 



((np) + (n a )) 



(3) 



= <(Am a ) 2 ), 

= (Am s Am i )/(r? s 77 i ). 



If QDEs rj s and r\i were known, Eq. (3) would give five 
constraints for the determination of six parameters M a 



and b a . This represents a serious problem arising from 
two points: (1) Only the first and second photoelectron 
moments are experimentally available with sufficient pre- 
cision and (2) The probe photon field is non-classical due 
to its predominantly paired character that enforces the 
use of at least six independent parameters in its realistic 
description. The solution of Eq. (3) can be expressed as 
a one-parameter system conveniently parameterized by 
the mean photon-pair number (n p ) (a = s,i): 

(n a ) = (m a )/ri a - (n p ), 



((Ah„) 2 ) = 
((An p ) 2 > = 



{(Am a ) 2 ) (Am s Am,) 1 - n a 



Va VsVi 

(Am s Ami)/(r) s r)i). 



2 \ m a), 



(4) 



However, also QDEs rj s and r\i are not known. As the 
third and higher photoelectron moments cannot be re- 
liably used for the chosen number of measurement rep- 
etitions [16], we suggest to minimize the declinations of 
experimental (/) and theoretical (p c ) photocount dis- 
tributions. A JPCD p c can be derived from a JPND p 
provided that the detection process is characterized. The 
JPND p of a field composed of pair, signal, and idler 
components can naturally be written as a convolution of 
three Mandel-Rice distributions [15]: 



min[n s .nt 

p(n s ,Ui) = ^2 p(n s - n;M s ,b s ) 

n=0 

x p(m - n; M h bi)p(n; M p , b p ); (5) 



p(n; M, b) = T(n + M)/[n! r(A/)]6"/(l + b) n+M . Follow- 
ing Eq. (2), numbers M a of modes and mean photon- 
numbers b a per mode can be derived from the first and 
second photon moments occurring in Eq. (4) as follows 

(a = p,s,i): 



((An a ) 2 ) 
(n a ) l! 



M a 



{(An a ) 2 ) - (n a )' 



(0) 



A detailed theory appropriate for a detector with AT 
pixels, QDE r\ and dark-count rate D = (d) /N has been 
developed in [16] and provided the probabilities T(c, n) 
of having c photocounts out of a field with n photons: 



T a (c,n) 



\ ^ / c 



a (i-A,)"-(i -va) n (-iy 



(-iy 



1=0 v 7 ^ ; 
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N a 1 - rj a 



,a = s, i. (7) 



Using the formulas in Eqs. (5) and (7) the JPCD p c can 
be expressed as: 

00 

Pc(c s ,Ci)= ^2 T s (c s ,n s )Ti(ci,ni)p(n s ,ni). (8) 

n s ,rii— 

Finally, the appropriate values of QDEs n s and rji and 
mean photon-pair number (n p ) are obtained by minimiz- 
ing the function T> quantifying the declinations of JPCD 
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Fig. 1. Minimum of declinations T> determined over the 
allowed values of (n p ) depending on QDEs rj s and n. 
In black areas close to rj s and n axes, Eqs. (3) have no 
solution. For optimum values of QDEs r} a = 0.243 and 
?]i = 0.235, JPND p(n s ,rii) is shown. 
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[pc{Cs,Ci) - f(c s ,d)Y 



(9) 



To practically demonstrate the power of the method, 
we discuss a typical measurement described in Fig. 1. 
Active area of the used iCCD camera is composed of 
128x128 independent pixels with equal QDEs obtained 
by hardware binning of the original 1 megapixel reso- 
lution. Processing the raw data eliminates spatial blur- 
ring of the detection spots. Signal and idler photons are 
captured in different areas of the photocathode. Some 
pixels in the photocathode are also reserved for moni- 
toring the noise. However, not all noises can be quanti- 
fied this way. For example, single photons arising from 
fluorescence in the nonlinear crystal are difficult to mon- 
itor. There also occur pump-intensity fluctuations that 
modify the photon-pair statistics. As an example, we 
analyze the measurement that has resulted in the fol- 
lowing photoelectron moments after eliminating dark 
counts: (m s ) = 2.411 ± 0.002, (m*) = 2.353 ± 0.004, 
((Am s ) 2 ) = 2.489 ± 0.003, {(Am;) 2 ) = 2.449 ± 0.005, 
and (Am s Ami) = 0.597±0.003. Applying the usual sim- 
plified method for the determination of QDE [neglecting 
noises in (3) and assuming ((An p ) 2 ) — (n p )], we arrive 
at the values r) s = (Am s Ami) / '( m i) = 0.254 ±0.001 and 
j]i = (Am s Am,)/(m s ) = 0.248 ±0.001. For comparison, 
covariance of the measured signal and idler photocounts 
was 0.238 ±0.001 [5]. 

On the other hand and analyzing the data along the 
developed method, we obtain the graph of minimum de- 
clinations T> depending on QDEs rj s and r\i as shown in 
Fig. 1. The plotted minimum values of declinations T> arc 
found after scanning over the mean photon-pair number 
(rip). The minimum value of declinations T> has occured 
for QDEs T} s = 0.243 ± 0.001 and r?; = 0.235 ± 0.001. 
Whereas efficiency r] s gives directly QDE of the camera, 
efficiency r/, encompasses also reflectivity of the mirror 
placed in the path of the idler beam in perfect agreement 
with its independently measured value. The comparison 
with the values written above shows that partitioning 
the noisy parts of the probe field results in lower values 



of QDEs. The method has also provided the complete 
characteristics of the probe field with 1% relative er- 
ror: b p = 0.058, M p = 170, b s = 33.2, M a = 0.0007, 
h = 10.6, and M* = 0.0101. Thus, the probe field has 
contained on average 9.91 pairs, 0.02 signal noise photons 
and 0.11 idler noise photons. Despite the fact that the 
probe field has been nearly-exclusivcly composed of pho- 
ton pairs (for the JPND, see Fig. 1), the effect of noise 
on the values of QDEs cannot be neglected. It holds that 
the larger the mean photon-pair number (n p ) the closer 
the values of QDEs obtained with the standard and the 
developed methods. 

In conclusion, we have developed a method for pre- 
cise determination of absolute detector efficiency of any 
photon-number resolving detector. The method allows 
to partition noise from the probe predominantly paired 
field thus providing, in principle, the precision in detector 
calibration limited only by the number of measurement 
repetitions. The improved precision and the possibility to 
use noisier probe fields make the method superior above 
other methods developed so far. 
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